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ABSTRACT:   This  paper  is  concerned  with  testing  normality  of  observations 
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1.   INTRODUCTION 

Numerous  tests  are  available  for  testing  the  normality  of  observa- 
tions (see  Mardia,  1980),  and  many  of  those  tests  can  be  modified  to 
test  the  normality  of  regression  disturbances  by  using  regression  resi- 
duals (see  Pierce  and  Kopecky,  1979).   Most  of  these  tests  have  been 
constructed  by  exploiting  some  unique  properties  of  the  normal  distri- 
bution that  are  not  shared  by  other  distributions.   It  is  not  known 
whether  they  possess  any  "optimal"  properties.   Recently,  3era  and 
Jarque  (1981),  3era  (1982),  Bera  and  John  (1983),  and  Poirier  et  al . 
(1984)  have  suggested  tests  employing  the  Lagrange  multiplier  (LM)  or 
Rao's  (1948)  score  principle  when  the  alternatives  belong  to  certain 
families  of  distributions  which  reduce  to  the  normal  distribution  under 
certain  parametric  restrictions.   Since  these  tests  are  based  on  the  LM 
principle  they  can  be  expected  to  possess  good  power  properties.   Simu- 
lation results  reported  in  Bera  and  Jarque  (1981)  reveal  that  the  test 
based  on  the  Pearson  alternatives,  henceforth  denoted  by  LM^,,  performs 
as  well  as  the  available  tests  in  finite  samples.   Statistic  LM^,  is 
very  easy  to  compute  (some  of  the  tests  which  have  been  found  to  have 
good  power  properties  are  not  very  easy  to  apply,  e.g.,  for  the  Shapiro 
and  Wilk  (1965)  test  we  require  certain  coefficients  which  are  not 
available  for  all  sample  sizes). 

One  drawback  of  the  Pearson  family  is  that  it  does  not  permit  the 
distributions  to  have  infinite  variance.   We  try  to  tackle  this  problem 
by  considering  the  stable  family  of  distributions.   Various  investiga- 
tions (see  for  instance,  Mandelbrot,  1963  and  Carlson,  1975)  have  shown 
that  certain  economic  variables  such  as  stock  price  changes  and  price 
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expectations  can  be  better  described  by  stable  distributions.   Another 
notable  aspect  of  this  family  is  that  the  distribution  of  the  regression 
disturbance  (which  is  assumed  to  be  the  sum  of  a  large  number  of  inde- 
pendent and  identically  distributed  random  variables  whose  individual 
effects  are  very  small)  may  belong  to  this  family  as  argued  by  Bartels 
(1977),  Koenker  and  Bassett  (1978)  and  others.   This  follows  from  the 
generalized  central  limit  theorem  which  states  that  for  a  distribution 
function  to  be  a  limit  distribution  of  a  sum  of  independent  and  identi- 
cally distributed  random  variables,  it  is  necessary  and  sufficient  that 
it  be  stable  (see  Gnedenko  and  Kolmogorov,  1954,  p.  162).   In  Section  2, 
we  attempt  to  derive  a  LM  test  for  normality  by  considering  the  stable 
family  as  the  alternative.   It  seems  that  no  computationally  simple 
test  can  be  developed  by  following  this  approach.   Therefore,  in  Section 
3,  we  investigate  whether  the  LM^  and  other  available  tests  can  detect 
non-normality  when  the  alternative  belongs  to  the  stable  family.   This 
is  done  through  a  simulation  study  which  closely  follows  the  design 
described  in  Bera  and  Jarque  (1981).   We  also  provide  a  table  from  which 
critical  points  of  USL,  for  different  sample  sizes  and  significance  levels 
can  easily  be  calculated.   Simulations  in  the  paper  were  carried  out  on 
the  Australian  National  University's  UNIVAC  1100  Computer  (Tables  I,  II 
and  III),  and  Cyber  175  (Table  IV)  at  the  University  of  Illinois.   The 
paper  is  concluded  in  Section  4  with  some  remarks. 

2.   LM  TEST  WITH  STABLE  ALTERNATIVES 

Let  us  consider  having  a  set  of  n  independent  observations 
u.,u-,...,u  ,  on  a  random  variable  u.   Assume  the  probability  density 
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function  (PDF)  of  u  is  a  member  of  the  stable  family.   Therefore,  we 
can  write  (see  Gnedenko  and  Kolmogorov,  1954,  p.  164) 

00 

f(u.)    =-^/      exp(itu    )    .    exp[iUt-y|t|a(l    +  i|j  u( t ,o) } Jdt        (1) 

J  —CO  "*  II 

—co  <  U  .   <C  °° 

j  ■  1,2, ... ,n 
where  i  =  /— 1,  t  is  any  real  number  and 

u(t,a)  =  tan  [^)    ,         if  a  *  1 

=  -  In 1 1  j  ,  if  a  =  1  . 

There  are  four  parameters:   the  location  parameter  y  e  (-ao,ao),  the  scale 
parameter  y  e  (0,°°),  the  skewness  parameter  6  e  [-1,1]  and  the  charac- 
teristic exponent  or  the  index  parameter  a  £  ( 0 , 2 J  .   The  stable  family 
includes  some  important  distributions  as  special  cases,  for  example,  the 

Cauchy  distribution  when  a  =  1  and  6=0;  the  normal  distribution  with 

2  2 

mean  p  and  variance  a      when  a  =  2,  6  =  0  and  y   =  a  /2.   The  distribution 

has  finite  variance  only  when  ct  =  2,  otherwise  the  variance  is  infinite. 

Since  we  are  interested  only  in  the  shape  of  the  distribution, 

without  loss  of  generality  we  can  assume  y  =  0  and  y  =  —.      Suppose  we 

also  assume  that  the  population  mean  exists,  i.e.,  a  >  1  then  we  can 

write  (1)  as 

00 

f(u  )   =     *    f     exp(itu  )   .    exp[-i   |t|a{i   +1"   tan(ZO)}]dt.      (2) 

J     2tt  j  2  t       2 » J 

—00 

A  test  for  normality  can  be  achieved  by  testing  H  :   a  =  2,  5  =  0  in 
(2).   However,  there  are  a  few  potential  difficulties  in  doing  so. 
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Under  H  ,  the  value  of  a  lies  on  the  boundary  of  the  parameter  space, 
and  in  addition,  at  a  =  2,  the  parameter  6  disappears  from  (2).   For 
simplicity  let  us  consider  the  symmetric  stable  distribution  by 
assuming  6=0.   Then  f(u.)  reduces  to 


f(u  )  -^J   cos(u.t)  exp[--^-]dt 


We  can  obtain  a  test  for  normality  by  testing  H  :   a  =  2  vs_  a  <  2,  based 
on  the  maximum  likelihood  estimate  (MLE)  of  a,  say  a.   This  is  suggested 
in  DuMouchei  (1971)  and  performance  of  this  test  has  been  investigated 
by  Saniga  et  al .  (1975).   However,  in  a  recent  paper,  DuMouchei  (1981) 
points  out  two  drawbacks  of  this  procedure.   First,  for  a  <  2,  a  is 
asymptotically  normal  with  mean  a  and  variance  equal  to  the  inverse  of 
the  Fisher's  information;  but  this  asymptotic  theory  fails  when  H  : 
a  =  2  is  true.   Second,  the  test  based  on  a  is  not  robust  to  the  assump- 
tion of  stability,  e.g.,  the  test  may  diagnose  the  exponential  distribu- 
tion as  having  infinite  variance.   (Of  course,  most  of  the  tests  will 
have  this  defect.)   Moreover,  the  test  suffers  from  computational  com- 
plexity. 

Let  us  now  explore  the  possibility  of  constructing  a  test  using  the 
Rao's  score  test  principle.   The  test  may  be  based  on  the  score  value 
3£(a)/3a  evaluated  at  a  =  2  where  £(a)  =  £._,  ln(f(u.)J.   More  specif i- 
cally,  we  examine  whether  the  following  quantity 


3Z(q)  m   3A(q) 
3a      3a 


n  T(u  ) 

=  e  — i- 

a=2  j-1  f<uj) 


where 
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3f(u .) 

f'(u.)  =     J 

J      9a 


a=2 


1  t   2 

=  -  —  /   cos(u.t)  exp  ( j)t   ln|t|dt 


and 


~  1 

f(u  )  =  -==  exp[ *■]  , 

J      i  It 


is  significantly  different  from  zero  using  a  one-sided  normal  test. 
By  numerical  integration,  3fc(a)/3a  and  an  estimate  of  its  variance 


(under  H  ) 


32l(q) 
3a2 


32£(a) 
3a2 


n  T'(u  )t(u.)-[?'(u.)J 

=  _  y.     - 1 = J 


a=2 


j-l 


[f(u.)J 


where 


3  f(u.) 

3a 


«  2.2 

38  ~  47  /   cos(u  t)exp(-  -^■)(ln|t|)2t2(-  -^  +  l)dt, 
a=2         -<» 


can  be  calculated.   However,  for  practical  purposes  the  calculations 
seem  to  be  quite  tedious.   One  reason  for  all  these  difficulties  is  that 
no  simple  expression  is  available  for  the  density  function  of  the  stable 
distribution. 

It  might  be  easier  to  construct  a  test  using  its  characteristic 
function  (CF)  which  has  a  simpler  expression.   For  example,  Epps  and 
Pulley  (1982)  suggested  a  test  based  on  a  weighted  integral  of  <j>n(t)  - 
^^(t)  where  $  (t)  and  4>n(t)  are  respectively  the  empirical  GF  and  CF  of 
a  normal  distribution.   To  obtain  high  power  against  the  stable  alterna- 
tives we  can  choose  the  value  of  t  to  maximize  the  power  instead  of 
taking  a  weighted  integral  over  t.   But  this  procedure  also  involves 
considerable  amount  of  computation. 
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3.   SIMULATION  STUDY 

Given  the  above  complexities  of  the  stable  distribution,  we  follow 
an  ad  hoc  approach.   In  particular,  we  examine  the  performances  of  the 
conventional  tests  against  the  stable  alternatives.   Many  of  these  tests 
are  based  on  the  second,  third  and  fourth  sample  moments,  but  for  the 
stable  distribution  the  population  counterparts  of  these  moments  do  not 
exist.   Therefore,  we  cannot  be  sure  about  the  distribution  of  these 
statistics  when  the  stable  alternative  is  true.   However,  if  we  find 
that  these  tests  can  detect  non-normality  when  the  alternatives  come 
from  stable  family,  for  practical  purposes  we  can  use  them  irrespective 
of  their  analytical  properties.   This  is  what  we  set  out  to  investigate 
in  this  section. 

Many  simulation  studies  have  been  done  on  assessing  the  performances 
of  the  conventional  tests  for  normality  against  stable  alternatives. 
For  instance,  Fama  and  Roll  (1971),  Saniga  et  al .  (1975),  Smith  (1975) 
and  Saniga  and  Hayya  (1977)  have  examined  the  powers  of  various  tests 
for  the  normality  of  observations  against  the  symmetric  stable  distri- 
butions, and  recently  Saniga  and  Miles  (1979)  have  extended  those 
studies  to  the  asymmetric  stable  family.   To  our  knowledge  there  is  no 
similar  systematic  investigation  on  testing  for  the  normality  of 
regression  disturbances.   Therefore,  for  our  simulation  experiment,  we 
consider  both  the  symmetric  and  asymmetric  stable  family  and  study  the 
powers  of  some  of  the  available  tests  for  the  normality  of  both  obser- 
vations and  regression  disturbances. 

We  consider  the  following  tests:   /b   b«,  D'Agostino  (1971)  D*, 
Pearson  et  al .  (1977)  R,  Shapiro  and  Wilk.  (1965)  W,  Shapiro  and  Francia 
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(1972)  Wf  and  Bera  and  Jarque  (1981)  LMM.   White  and  MacDonald  (1980) 

IN 

describe  the  first  six  tests  and  LM^,  takes  the  following  form 

(/b  )2   (b  -3)2 
U\   =  n[—6 +  W~]  "  (3) 

A  short  derivation  of  this  statistic  is  given  in  the  Appendix;  for 
details  see  Bera  and  Jarque  (1981).   Under  normality  (H  ),  LM  is 
asymptotically  distributed  as  a  Xo*   Distributions  under  the  alternative 
hypothesis  are  generated  from  the  stable  family  using  the  density  func- 
tion given  in  (2)  with  a  =  1.0(.3)1.9  and  5  -  0.0(. 25)1.0.   To  do  this 
we  make  use  of  the  computer  program  described  in  Chambers  et  al .  (1976). 
Observations  under  normality  (a  =  2.0,  <5  =  0.0)  are  generated  using  the 
IBM  subroutines  RANDO  and  RANDU. 

To  calculate  powers  for  n  =  20 ,  35,  50  and  100  and  tests  /b, ,  b~ » 
D*  and  R,  we  use  the  critical  points  given  in  White  and  MacDonald  (1980, 
p.  20).   For  n  =  200,  significance  points  for  /b,  ,  b„  and  D*  are  obtained 
respectively  from  Pearson  and  Hartley  (1962,  p.  183),  D'Agostino  and 
Pearson  (1973,  p.  615)  and  D'Agostino  (1971,  p.  343);  and  for  the  R  test 
we  extrapolated  the  points.   For  W,  W'  and  LMV,  we  compute  the  signifi- 
cance points  by  simulation  using  500  replications,  in  order  to  obtain  an 
empirical  significance  level  equal  to  .10.   For  example,  for  a  given  n, 
the  significance  point  for  W  is  W(50)  -  the  50th  largest  of  the  values 
of  W  in  the  500  replications  under  normal  observations  and  similarly  for 
W'.   For  LM„,  the  point  is  set  as  LM^(450).   With  the  empirical  signi- 
ficance level  set  at  .10,  easier  power  comparisons  among  the  one-sided 
tests  W,  W'  and  UL.   can  be  made.   Since  /b  ,  b  and  D*  are  two-sided 
tests  and  R  is  a  four-sided  test,  it  is  difficult  to  adjust  the  signi- 
ficance points  to  achieve  an  empirical  significance  level  equal  to  .10. 
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Estimated  powers  of  each  test  (obtained  by  dividing  the  number  of 
times  H  is  rejected  by  500)  for  testing  for  the  normality  of  observa- 
tions are  reported  in  Table  I  for  n  =  20,  35,  50,  100  and  200,  except 
for  W  which  cannot  be  computed  for  n  >  50.   The  maximum  standard  error 


of  each  entry  in  the  table  is  / .5(l-.5)/500  -    .022.   In  the  table,  the 
highest  power  is  underlined  (for  each  alternative  and  sample  size) 
except  when  three  or  more  tests  have  this  power.   To  save  space,  results 
corresponding  to  only  nine  combinations  of  a  and  6  namely,  a  =  1.0,  1.6 
and  1.9  and  5  =  0.0,  0.5  and  1.0  are  given.   They  are  representative  of 
the  overall  results. 

For  each  n,  the  first  row  gives  the  estimated  sizes  of  the  tests: 
/b,  ,  b?,  D*  and  R;  they  are,  in  most  cases,  around  .10.   Therefore, 
power  comparison  of  these  tests  with  W,  W,  and  LM„  is  valid.   From  the 
table  we  note  that  for  a  given  value  of  a  the  powers  are  not  sensitive 
to  changes  in  the  value  of  the  skewness  parameter  <5 .   However,  when  5 
is  kept  fixed  and  a  is  varied,  powers  change  considerably.   The  overall 
behavior  of  the  first  five  tests,  /b  ,  b  ,  D* ,  R  and  W  is  very  similar 
to  those  reported  in  Saniga  and  Miles  (1979).   Two  additional  tests  we 
consider  are  W  and  LM^,,  and  it  is  seen  that  they  outperform  the  remain- 
ing ones.   In  particular,  LM  does  very  well  and  it  has  the  highest 
power  in  33  cases  out  of  45  cases  considered;  when  it  does  not,  its 
power  is  within  less  than  2  percent  of  the  maximum  power.   It  has  very 
good  power  for  smaller  sample  sizes  and  also  when  a  is  near  to  2,  i.e., 
when  the  alternative  is  close  to  the  normal  distribution. 

We  now  study  the  power  of  the  tests  for  the  normality  of  unobserved 
regression  disturbances.   The  tests  we  consider  are  the  same  as  discussed 


20 


35 


50 


100 


200 
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TABLE    I 


Estimated    powers    of    the    tests    against   stable   alternatives 


/b2  b2  D*  R  W  W*  LMjj 


2.0 

0.0 

.076 

.062 

.096 

.080 

.100 

.100 

.100 

1.0 

0.0 

.804 

.872 

.902 

.872 

.896 

.912 

.898 

1.0 

0.5 

.832 

.842 

.890 

.886 

.902 

.922 

.910 

1.0 

1.0 

.956 

.806 

.894 

.938 

.968 

.966 

.968 

1.6 

0.0 

.486 

.462 

.466 

.482 

.472 

.518 

.538 

1.6 

0.5 

.504 

.466 

.476 

.500 

.496 

.536 

.542 

1.6 

1.0 

.558 

.422 

.436 

.518 

.536 

.550 

.578 

1.9 

0.0 

.234 

.196 

.192 

.232 

.204 

.222 

.244 

1.9 

0.5 

.220 

.190 

.204 

.216 

.206 

.222 

.224 

1.9 

1.0 

.212 

.172 

.188 

.190 

.204 

.222 

.246 

2.0 

0.0 

.082 

.100 

.090 

.086 

.100 

.100 

.100 

1.0 

0.0 

.880 

.982 

.996 

.982 

.982 

.994 

.986 

1.0 

0.5 

.908 

.972 

.990 

.970 

.974 

.994 

.982 

1.0 

1.0 

.992 

.960 

.982 

.994 

.994 

.992 

.992 

i.6 

0.0 

.618 

.692 

.698 

.684 

.636 

.720 

.732 

1.6 

0.5 

.672 

.672 

.688 

.690 

.634 

.740 

.730 

1.6 

1.0 

.810 

.660 

.682 

.786 

.776 

.788 

.812 

1.9 

0.0 

.298 

.268 

.266 

.290 

.230 

.310 

.312 

1.9 

0.5 

.284 

.272 

.266 

.290 

.250 

.314 

.308 

1.9 

1.0 

.286 

.256 

.260 

.292 

.292 

.306 

.306 

2.0 

0.0 

.074 

.094 

.074 

.086 

.100 

.100 

.100 

1.0 

0.0 

.930 

.988 

.992 

.990 

.986 

.994 

.992 

1.0 

0.5 

.938 

.988 

.994 

.994 

.992 

.996 

.996 

1.0 

1.0 

1.000 

.986 

.998 

1.000 

1.000 

1.000 

1.000 

1.6 

0.0 

.712 

.806 

.808 

.802 

.692 

.828 

.850 

1.6 

0.5 

.740 

.798 

.804 

.810 

.714 

.830 

.836 

1.6 

1.0 

.856 

.716 

.742 

.834 

.818 

'  .840 

.860 

1.9 

0.0 

.304 

.316 

.306 

.322 

.224 

.346 

.354 

1.9 

0.5 

.330 

.328 

.304 

.328 

.256 

.354 

.374 

1.9 

1.0 

.386 

.322 

.316 

.368 

.308 

.372 

.406 

2.0 

0.0 

.068 

.108 

.114 

.112 

.100 

.100 

1.0 

0.0 

.946 

1.000 

1.000 

1.000 

1.000 

1.000 

1.0 

0.5 

.980 

1.000 

1.000 

1.000 

1.000 

1.000 

1.0 

1.0 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.6 

0.0 

.828 

.954 

.950 

.938 

.952 

.970 

1.6 

0.5 

.870 

.948 

.952 

.950 

.968 

.966 

1.6 

1.0 

.992 

.938 

.952 

.990 

.992 

.992 

1.9 

0.0 

.416 

.462 

.452 

.470 

.498 

.514 

1.9 

0.5 

.420 

.460 

.438 

.468 

.510 

.500 

1.9 

1.0 

.528 

.444 

.426 

.498 

.520 

.556 

2.0 

0.0 

.102 

.106 

.108 

.110 

.100 

.100 

1.0 

0.0 

.984 

1.000 

1.000 

1.000 

1.000 

1.000 

1.0 

0.5 

.980 

1.000 

1.000 

1.000 

1.000 

1.000 

1.0 

1.0 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.6 

0.0 

.892 

.998 

.998 

.996 

.996 

.996 

1.6 

0.5 

.932 

.99~6 

.99"3" 

.994 

.998 

.998 

1.6 

1.0 

1.000 

.994 

.996 

1.000 

1.000 

1.000 

1.9 

0.0 

.522 

.652 

.632 

.668 

.696 

.696 

1.9 

0.5 

.620 

.664 

.624 

.672 

.714 

.710 

1.9 

1.0 

.786 

.694 

.682 

.774 

.774 

.800 
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above  but  we  compute  them  using  the  ordinary  least  squares  (OLS)  resi- 
duals rather  than  the  true  disturbances  u. .   The  result  of  Pierce  and 
Kopecky  (1979)  indicates  that  this  modification  is  valid  asymptotically. 
We  denote  the  modified  statistics  by  /b,,  b„  etc.  and  consider  the  linear 
regression  model 


y .  =  I   S.x..  +  u.  j  =  1,2, ...,n. 

3        iasl  i  Ji    J  J 

Simulation  results  studying  the  relative  powers  of  tests  for  normality 

of  u  =  (u, ,U2,...,u  )',  computed  using  the  OLS  residuals  u  = 

(u  ,u?,...,u  )',  depend  on  the  design  matrix  X  since  u  =  Mu  where  M  is 

an  n  x  n  matrix  defined  by  M  =  I  -  X(X'X)   X1.   Therefore,  we  consider 

two  sets  of  regressors.   The  first  set  is  defined  as  in  White  and 

MacDonald  (1980,  p.  20),  i.e.,  we  set  x   =  1  (j  =  1,2,... ,n),  and 

generate  X~  ,  X«  and  X,  from  a  uniform  (U)  distribution.   The  last  three 

regressions  are  transformed  to  have  mean  zero  and  variance  25.   The 

specific  values  of  the  means  and  variances  of  these  regressors  have  no 

effect  on  the  simulation  results.   For  the  second  set,  we  set  x.,  =1 

jl 

(j  =  l,2,...,n),  and  generate  X2  from  a  N(10,25),  X-  from  a  U[7.5,  12. 5] 

2 
and  X,  from  a  y  ,  ,^  • 
4        A10 

The  results  corresponding  to  these  two  data  sets  are  given  in  Tables 
II  and  III.   To  calculate  the  powers  we  utilize  the  same  significance 
points  used  in  the  calculation  of  powers  in  Table  I,  except  for  W,  W' 
and  LM  ,  for  which  we  use  empirical  10  percent  critical  points.   For 

si 

instance,  the  significance  point  for  W  is  W(50). 

The  first  point  to  note  from  these  tables  is  that  in  almost  all  cases 
modified  tests  give  lower  powers  than  the  corresponding  tests  using 
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TABLE  II 


Estimated  powers  of  the  tests  against  stable  alternatives 
using  the  OLS  residuals 


(Regressors:      X. 

-  i,x2 

>VX4  - 

U) 

a 

6 

*x 

~b2 

D* 

R 

w 

W' 

^ 

2.0 

0.0 

.080 

.090 

.112 

.092 

.100 

.100 

.100 

1.0 

0.0 

.734 

.778 

.792 

.794 

.782 

.822 

.816 

1.0 

0.5 

.760 

.754 

.772 

.788 

.768 

.814 

.830 

1.0 

1.0 

.876 

.756 

.786 

.860 

.866 

.874 

.884 

■  ■  2°    1:1 

0.0 
0.5 

.402 
.408 

.348 
.366 

.362 
.352 

.368 
.398 

.370 
.396 

.414 
.446 

.434 
.448 

1.6 

1.0 

.486 

.366 

.384 

.450 

.456 

.474 

.508 

1.9 

0.0 

.198 

.186 

.190 

.184 

.188 

.208 

.200 

1.9 

0.5 

.200 

.188 

.184 

.184 

.180 

.192 

.208 

1.9 

1.0 

.174 

.192 

.176 

.192 

.184 

.172 

.194 

2.0 

0.0 

.096 

.098 

.102 

.082 

.100 

.100 

.100 

1.0 

0.0 

.854 

.966 

.976 

.964 

.948 

.974 

.968 

1.0 

0.5 

.904 

.966 

.976 

.970 

.960 

.982 

.980 

1.0 

1.0 

.988 

.940 

.964 

.984 

.992 

.992 

.988 

1c          1.6 
n  =  35          1.6 

0.0 
0.5 

.612 
.634 

.632 
.614 

.646 
.638 

.630 
.654 

.582 
.600 

.672 
.684 

.696 
.686 

1.6 

1.0 

.758 

.594 

.614 

.722 

.696 

.742 

•  .754 

1.9 

0.0 

.258 

.236 

.248 

.252 

.238 

.276 

.282 

1.9 

0.5 

.266 

.258 

.244 

.270 

.248 

.286 

.296 

1.9 

1.0 

.268 

.226 

.230 

.266 

.258 

.280 

.276 

2.0 

0.0 

.086 

.090 

.074 

.066 

.100 

.100 

.100 

1.0 

0.0 

.906 

.982 

.990 

.982 

.978 

.984 

.984 

1.0 

0.5 

.932 

.986 

.992 

.994 

.988 

.996 

.994 

1.0 

1.0 

1.000 

.980 

.992 

.998 

.998 

.998 

.998 

n  =  50          1.6 

0.0 

.688 

.772 

.786 

.780 

.682 

.794 

.808 

0.5 

.744 

.766 

.782 

.802 

.708 

.800 

.814 

1.6 

1.0 

.842 

.696 

.708 

.810 

.774 

.824 

.838 

1.9 

0.0 

.290 

.310 

.294 

.302 

.232 

.308 

.324 

1.9 

0.5 

.322 

.330 

.314 

.326 

.236 

.328 

.352 

1.9 

1.0 

.348 

.312 

.290 

.338 

.302 

.366 

.388 

2.0 

0.0 

.066 

.098 

.098 

.090 

.100 

.100 

1.0 

0.0 

.936 

1.000 

1.000 

1.000 

1.000 

1.000 

1.0 

0.5 

.978 

1.000 

1.000 

1.000 

1.000 

1.000 

1.0 

1.0 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

n  -  ioo     };* 

1.6 

0.0 

.818 

.948 

.940 

.934 

.946 

.960 

0.5 
1.0 

.866 
.986 

.936 
.926 

.942 
.934 

.940 
.986 

.958 
.986 

.958" 
.986 

1.9 

0.0 

.410 

.452 

.430 

.458 

.484 

.496 

1.9 

0.5 

.432 

.444 

.432 

.470 

.486 

.500 

1.9 

1.0 

.520 

.436 

.418 

.496 

.528 

.548 

2.0 

0.0 

.092 

.122 

.114 

.100 

.100 

.100 

1.0 

0.0 

.980 

1.000 

1.000 

1.000 

1.000 

1.000 

1.0 

0.5 

.982 

1.000 

1.000 

1.000 

1.000 

1.000 

1.0 

1.0 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

n-200       J;J 

0.0 

.888 

.998 

.998 

.998 

.996 

.996 

0.5 

.924 

.996 

.996 

.994 

.996 

.996 

1.6 

1.0 

1.000 

.994 

.994 

1.000 

1.000 

1.000 

1.9 

0.0 

.512 

.652 

.638 

.660 

.674 

.678 

1.9 

1.9 

0.5 

1.0 

.620 
.772 

.650 
.686 

.630 
.682 

.664 

.770 

.694 
.778 

.694 
.800 
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TABLE  III 


estimated  powers  of  the  tests  against  stable  alternatives 
using  the  OLS  residuals 

(Regressors:   Xj  -  1,  Xg  ~  M,  Xj  ~  U,  X^  ~  X^Q) 


a 

6 

*1 

~b2 

D* 

• 

s 

W 

W' 

^ 

2.0 

0.0 

.080 

.104 

.106 

.088 

.100 

.100 

.100 

1.0 

0.0 

.726 

.766 

.780 

.782 

.764 

.802 

.822 

1.0 

0.5 

.766 

.740 

.778 

.790 

.770 

.812 

.822 

1.0 

1.0 

.864 

.700 

.746 

.826 

.854 

.862 

.868 

» ■ «    1.1 

0.0 

.396 

.404- 

.418 

.416 

.390 

.444 

.452 

0.5 

.432 

.396 

.404 

.426 

.414 

.446 

.484 

1.6 

1.0 

.450 

.338 

.370 

.414 

.442 

.456 

.480 

1.9 

0.0 

.202 

.162 

.170 

.184 

.176 

.170 

.206 

1.9 

0.5 

.202 

.166 

.158 

.186 

.178 

.182 

.206 

1.9 

1.0 

.196 

.176 

.176 

.194 

.226 

.242 

.234 

2.0 

0.0 

.068 

.092 

.084 

.088 

.100 

.100 

.100 

1.0 

0.0 

.870 

.958 

.974 

.960 

.952 

.972 

.964 

1.0 

0.5 

.900 

.960 

.976 

.966 

.966 

.982 

.970 

1.0 

1.0 

.988 

.932 

.966 

.984 

.986 

.990 

.986 

n   =  35          1.6 

0.0 

.602 

.648 

.646 

.658 

.584 

.676 

.694 

0.5 

.646 

.622 

.648 

.650 

.626 

.680 

."694" 

1.6 

1.0 

.776 

.602 

.630 

.734 

.720 

.750 

.776 

1.9 

0.0 

.296 

.260 

.256 

.262 

.258 

.310 

.306 

1.9 

0.5 

.294 

.268 

.258 

.290 

.266 

.310 

.316 

1.9 

1.0 

.284 

.236 

.244 

.264 

.290 

.304 

.312 

2.0 

0.0 

.062 

.074 

.080 

.066 

.100 

.100 

.100 

1.0 

0.0 

.910 

.982 

.990 

.982 

.976 

.986 

.986 

1.0 

0.5 

.930 

.988 

.992 

.992 

.990 

.996 

.996 

1.0 

1.0 

1.000 

.978 

.996 

.998 

.998 

.998 

1.000 

„.          1.6 
n   "   50          1.6 

0.0 

.684 

.758 

.780 

.766 

.678 

.798 

.810 

0.5 

.714 

.760 

.770 

.780 

.688 

.796 

.810 

1.6 

1.0 

.842 

.686 

.722 

.804 

.792 

.840 

.848 

1.9 

0.0 

.280 

.280 

.270 

.294 

.230 

.332 

.326 

1.9 

0.5 

.332 

.300 

.298 

.318 

.244 

.338 

.366 

1.9 

1.0 

.340 

.294 

.294 

.330 

.294 

.368 

.370 

2.0 

0.0 

.070 

.096 

.090 

.078 

.100 

.100 

1.0 

0.0 

.940 

1.000 

1.000 

1.000 

1.000 

1.000 

1.0 

0.5 

.976 

1.000 

1.000 

1.000 

1.000 

1.000 

1.0 

1.0 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

n-ioo     };66 

0.0 
0.5 

.810 
.870 

.948 
.942 

.944 
.944 

.928 
.948 

.948 
.966 

.962 
.964 

1.6 

1.0 

.986 

.936 

.938 

.986 

.986 

.988 

1.9 

0.0 

.416 

.464 

.436 

.466 

.478 

.510 

1.9 

0.5 

.428 

.458 

.432 

.472 

.498 

.512 

1.9 

1.0 

.516 

.438 

.408 

.498 

.514 

.552 

2.0 

0.0 

.068 

.112 

.130 

.124 

.100 

.100 

1.0 

0.0 

.986 

1.000 

1.000 

1.000 

1.000 

1.000 

1.0 

0.5 

.976 

1.000 

1.000 

1.000 

1.000 

1.000 

1.0 

1.0 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

n   -   200        }*S 

I  .0 

0.0 

.882 

.998 

.998 

.998 

.998 

.998 

0.5 

.926 

.996 

.996 

.994 

.998 

.998 

1.6 

1.0 

1.000 

.996 

.996 

1.000 

1.000 

1.000 

1.9 

0.0 

.520 

.648 

.632 

.662 

.694 

.688 

1.9 

0.5 

.622 

.664 

.638 

.670 

.708 

.710 

1.9 

1.0 

.774 

.678 

.666 

.766 

.774 

.790 
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the  original  disturbances  and  that  these  power  differences  diminish  as 
n  increases.   This  is  easily  seen  by  comparing  the  powers  in  Table  I 
with  those  in  Tables  II  and  III.   This  illustrates  the  observed  tendency 
of  the  OLS  residuals  to  be  "more  normal"  than  the  true  disturbances. 
Gnanadesikan  (1977,  p.  265)  calls  this  tendency  the  "supernormality"  of 
the  residuals.   Bassett  and  Koenker  (1982)  provide  a  simple  explanation 
of  this  behavior.   Each  residual  can  be  decomposed  into  true  disturbance 
and  a  term,  which  behaves  asymptotically  as  normal  irrespective  of  the 
true  nature  of  the  disturbance,  pushing  the  distribution  of  the  OLS  resi- 
duals towards  normality.   We  did  not  make  any  attempt  to  use  Theil's 
(1965)  best  linear  unbiased  scalar  (BLUS)  residuals  since  the  Monte 
Carlo  results  of  Huang  and  flolch  (1974)  show  the  superiority  of  OLS  over 
BLUS  residuals  for  testing  the  normality  of  the  regression  disturbances. 
We  also  find  that  the  ranking  of  the  tests  in  these  tables  is  approxi- 
mately the  same  as  in  Table  I.   LiyL.  has  the  highest  power  in  around  75 

percent  of  the  cases.   After  LM„,  the  next  best  test  is  W'. 

2 
Earlier  we  noted  that  LM^  is  asymptotically  distributed  as  a  Xo* 

However,  for  small  samples,  the  actual  critical  values  can  be  very 

2 
different  from  the  x?  critical  points.   So  the  use  of  asymptotic  criti- 
cal values  would  lead  to  misleading  inference  in  small  samples.   In 

Table  IV  we  report  the  ratios  of  empirical  mean  values  LM^  under  H  and 

2 
the  corresponding  Xo  mean  which  is  equal  to  2,  for  some  selected  values 

of  n.   For  this  table  observations  under  normality  are  generated  using 

the  IMSL  subroutine  GGNML.   The  empirical  mean  values  are  based  on 

10,000  replications.   Therefore,  the  maximum  standard  error  of  each 


entry  will  be  /1/10000  =  .01  and  we  can  expect  the  true  ratios  to  be 
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TABLE  IV 


Ratio  of  empirical  and  asymptotic  mean  values  of 
the  LM  test  statistic 


Sample  size 

Ratio 

10 

.46420 

20   . 

.65435 

30 

.73690 

40 

.77990 

50 

.83260 

60 

.84245 

70 

.86580 

80 

.89150 

Sample  size 

Ratio 

90 

.90105 

100 

.91255 

120 

.91560 

140 

.93845 

160 

.94680 

180 

.97440 

200 

.98145 

400 

.98620 
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within  +  .02  of  the  reported  values.   From  the  table  critical  values 

for  any  given  significance  level  and  sample  size  can  be  obtained  by 

2 
multiplying  the  given  ratio  with  the  x?  critical  value.   Let  n  =  30 

2 
and  the  significance  level  is  .01.   The  X2  critical  value  is  9.21,  and 

hence  it  would  be  more  appropriate  to  use  .73690  x  9.21  =  6.79  as  the 

critical  value  for  the  sample  size  30.   For  intermediate  values  of  n, 

the  ratios  can  be  obtained  by  interpolation. 


4.   CONCLUSION 

In  this  paper  we  have  investigated  the  performances  of  some  of  the 
available  tests  for  normality  against  stable  alternatives  by  simulation, 
It  has  been  observed  that  the  test  suggested  in  3era  and  Jarque  (1981) 
and  Shapiro  and  Francia  (1972)  W  test  are  most  effective  in  detecting 
non-normality  of  both  observations  and  regression  disturbances.   It 
seems  no  computationally  simple  test  can  be  developed  using  the  LM 
principle  by  directly  considering  the  density  function  of  the  stable 
family.   Possibly,  a  simple  test  can  be  constructed  from  its  charac- 
teristic function.   Further  work  in  this  direction  will  be  fruitful. 
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APPENDIX 

Let  £.(0)  denote  the  log-density  function  for  the  i-th  observation 
where  0  is  the  parameter  vector.   If  we  have  n  independent  observations, 
the  log-likelihood  function  is  Z(0)  =  I.  ,  £.(0).   Let  0  denote  the  MLE 
under  the  null  hypothesis  H  .   Then  the  Rao's  score  or  the  LM  test  sta- 
tistic to  test  H  is  given  by 

LM  =  d»(0)r1(0)d(0)  (A.l) 

where  d(8)  =  di(B)/dd,    1(9)  =  -E[3  £(e)/8030'J>  and  "~"  denotes  the 
expressions  are  evaluated  at  0  =  0 .   We  assume  that  the  PDF  of  u  is  a 
member  of  the  Pearson  family,  i.e., 

c, -u. 

exp[/ J du  ] 

c  -c  u  +c_u 


f(u.)  = °   V   2  1 

j     oo         c  — u 

/  exp[/ dujdu 


2 

C  -CtU+C0U 
o   1    2 


-co  <  U    K    °° 

J 

J  =  1 , 2 , • • • , n« 


When  c,  =  c~  =  0,  f(u.)  reduces  to  a  normal  density  function  with 
1    2         y 

mean  zero  and  variance  c  .   Therefore,  we  test  H  :   0  =  (c  ,c,,c0)' 

o  o         o  1   2 

(cQ,  0  0)'.   Let  us  denote 

c,-u 
h(9;u)  =  / du. 


co"clu+c2u 


Then 


Jt(0)  =  -n.ln[/  exp{h(0;u)}J  +  E  h(9;u  ). 
-»  j»l     J 

Using  this  expression  of  £(0)  we  can  show  that 
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u3  ul        M4  3 

d»<e)  =  q[o, r+TT'  — h  -v 

3u2Z        U2     4u2Z 


and 


1(0)    =   a 


2y 


4i 


3y 


2u, 


where  y.    =  E .,    u./n,    i   =  1,2,3,4   and  9    =    (p?,    0,    0)'.      Substituting 


d(9)    and    1(9)    in    (A.l),    we   obtain 

2 


3y 


M3  1      u4  2  Jyl  u3yl 

3        24  2  l\i  >  2 

6p2  P2  2  u2 

If    observations    are   measured    from   the   sample   mean,    y,    =  0,    and    then   LM 

2        3 
reduces    to    UL,  given   in   equation    (3)    by   defining   b,    =  y-   /u2     an<^   ^o   = 

/      2 
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